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We study diffusion of the Josephson phase in the asymmetric SQUID subjected to a time-periodic current and 
pierced by an external magnetic flux. We analyze a relation between phase diffusion and quality of transport 
characterized by the dc voltage across the SQUID and efficiency of the device. In doing so, we concentrate 
on the previously reported regime [J. Spiechowicz and J. Luczka, New J. Phys. 17, 023054 (2015)] for which 
efficiency of the SQUID attains a global maximum. For long times, the mean-square displacement of the phase 
is a linear function of time, meaning that diffusion is normal. Its coefficient is small indicating rather regular 
phase evolution. However, it can be magnified several times by tailoring experimentally accessible parameters 
like amplitudes of the ac current or external magnetic flux. Finally, we prove that in the deterministic limit 
this regime is essentially non-chaotic and possesses an unexpected simplicity of attractors. 
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A Superconducting quantum interference device 
(SQUID) is one of the most important element of 
apparatuses of research laboratories worldwide. 
It exhibits a wide variety of phenomena and has 
been successfully used for testing the fundamen¬ 
tals of quantum mechanics, quantum information 
and chaotic phenomena in classical physics. In 
the semiclassical regime, dynamics of the Joseph¬ 
son phase in the SQUID can be visualized as a 
Brownian particle moving in a periodic potential. 
This correspondence has allowed to study trans¬ 
port properties of the SQUID by applying meth¬ 
ods of Langevin equation. Earlier research in this 
field has been concentrated on the influence of 
thermal noise and external driving on current- 
voltage characteristics. The present paper stud¬ 
ies diffusion of the Josephson phase. We consider 
a special type of the SQUID which operates as a 
Brownian ratchet. In particular, we analyze the 
phase diffusion in the regime for which the effi¬ 
ciency of the SQUID is globally maximal. 


I. INTRODUCTION 

The SQUIDs 1 are elements of ultrasensitive electric 
and magnetic measurement systems. They are also ex¬ 
ploited in a number of commercial applications in in¬ 
dustrial metrology, geophysical systems and medicine 
as a noninvasive technique for investigating the human 
bodyPEO In terms of the Stewart-McCumber model,®^ 
both the dynamics of the Josephson phase in the SQUID 
and the position of a Brownian particle can be described 
by similar Langevin equation. 6 Therefore time evolution 


^Electronic mail: jerzy.luczka@us.edu.pl 


of the Josephson phase is analogous to a random mo¬ 
tion of a Brownian particle in a spatially periodic poten¬ 
tial. While prior studies have mostly been focused on the 
current-voltage characteristics, diffusion of the Josephson 
phase, i.e. its mean-square displacement, has not been in¬ 
tensively analyzed. There are only several papers closely 
or loosely related to this latter subjects 11 

In this work, we study diffusion of the Josephson 
phase in the asymmetric squieGshui 

which is composed 

of three capacitively and resistively shunted Josephson 
junctions. Two junctions are collocated in series in one 
half-piece of the ring and the third junction is disposed 
in the other half of the ring. The SQUID is driven by a 
time-periodic current and subjected to an external con¬ 
stant magnetic flux. Its dynamics is extremely rich and 
complex even in the deterministic case, involving har¬ 
monic, subharmonic, quasiperiodic and chaotic trajec¬ 
tories. At non-zero temperature, thermal fluctuations 
lead to diffusive behaviour with random escape events 
among possibly coexisting attractors and with the sys¬ 
tem typically exploring its whole phase space. Due to 
a multi-dimensionality of the parameter space describing 
the model it would be extremely difficult to perform a 
complete analysis of diffusion process of the Josephson 
phase in such a setup. Therefore we consider it only in 
the previously reported regime 19 for which the Stokes ef¬ 
ficiency of the SQUID attains its global maximum. 

The paper is structured as follows. In Sec. [TT] we 
present a model of the asymmetric SQUID in terms of 
the Langevin equation for the Josephson phase. Next, 
we define several quantifiers characterizing the quality of 
transport process occurring in this setup. Among them 
are those that relate directly to an asymptotic long time 
stationary average voltage drop like its variance and the 
Stokes efficiency of the SQUID. However, we also study 
in this context those connected with the realization of 
stochastic motion of the Josephson phase like its mean 
squared displacement and a diffusion coefficient. We 
also combine these complementary quantities into the so 
called Peclet number. In Sec. |IV|we present results based 






FIG. 1. Schematic asymmetric SQUID composed of three 
Josephson junctions and the equivalent circuit built from two 
junctions. The Josephson phase difference is tpi = + ifd, 

the externally applied current is I, the external magnetic 
flux is $ e and the instantaneous voltage across the SQUID 
is V = V(t). 


on the comprehensive numerical simulation of the stud¬ 
ied system in a regime for which the transport efficiency 
is globally maximal. In particular, we demonstrate a 
possibility of steering of diffusive behaviour of the phase 
by tuning the experimentally accessible parameters like 
the amplitude of ac current applied to the system or the 
external magnetic flux. It turns out that the diffusion 
coefficient can be significantly enhanced by small varia¬ 
tion of these quantities. Moreover, we establish a clear 
relation between the directed transport, its efficiency and 
the diffusive motion of the Josephson phase. It is very 
important that the regime of globally maximal transport 
quality measured by the Stokes efficiency is also effective 
in the sense that phase motion is ordered and regular. Fi¬ 
nally, the last section provides summary and conclusions. 


II. MODEL OF THE ASYMMETRIC SQUID 

The asymmetric SQUIEG®E3is presented in Fig. [l] It 
is formed by a superconducting loop with two resistively 
and capacitively shunted Josephson junctioiifP^l j n the 
left arm and only one in the right arm. The crosses de¬ 
note the junctions and tpk = (k = u , d , 1,2) are the 
phase differences across them. Each junction is charac¬ 
terized by the capacitance Ck, resistance Rk and critical 
Josephson current Jk , respectively. To reduce a number 
of parameters of the model, we consider only a special 
case when two junctions in the left arm are identical, 
i.e. J u = Jd = JuRu = Rd = Ri/2,C U — Cd = 2C \. 
In some regimes 2 -' they can be considered as one for 
which the supercurrent-phase relation takes the form 
Ji sin (ipi/2), where pi = (p u + ipd- Additionally, the 
SQUID is threaded by an external magnetic flux <b e . As 
a consequence the effective potential experienced by the 
phase ipi forms a ratchet structure.'- 2 The device is driven 
by an external current I = I(t) which is composed of the 
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FIG. 2. The potential Q for the symmetric case j = 0 (solid 
red line) in comparison with the ratchet potential for j = 
1/2 and two values of the external magnetic flux 4> e = 7 t/ 2 
(dashed green line) and $ e = —7 t/ 2 (dotted blue line). 

static dc current Iq and the ac component of amplitude 
A and angular frequency fl, namely, 

I(t) = Io + Acos(flt). (1) 

The Langevin equation for the phase ip = is of the 

forrrPU 

hi. , . , . /2fcsT 

te Ctp+ 2- e RV + J (.v) = m + ^j^ r m, ( 2 ) 

where the effective supercurrent J(<p) reads 

= J isin(!) + J 2 sin(iyj +$ e ). (3) 

The parameters are: C = C\ + C 2 , R~ l = Ri 1 + R^ 1 , 
ks is the Boltzmann constant, T is temperature of the 
system and >I> e = 27r<I> e /<i>o is the dimensionless external 
magnetic flux. Thermal fluctuations are modelled by 5- 
correlated Gaussian white noise £ (t) of the statistics: 

m) = 0 , mas)) = s(t~ S ). ( 4 ) 

The average voltage across the device can be calculated 
from the relation 

(V) = (h/2e)(ipi) = (h/2e)(<p 2 ), (5) 

where the averaging is over the period 27r/Q of the ac cur¬ 
rent. This relation is also valid in the long time regime 
when additionally averaging is performed over initial con¬ 
ditions and all realizations of thermal noise. 

It is useful to interpret Eq. ([2]) in the mechanical 
framework as a model of the inertial Brownian particle 
subjected to the conservative force J(tp) and propelled by 
the time dependent driving I(t). In this correspondence, 
the particle position x translates to the phase ip, its ve¬ 
locity v = x to the voltage V, the conservative force to 
the supercurrent the external force to the current 
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J(f), the mass m to the capacitance C and the friction 
coefficient 7 to the normal conductance G — 1/R. 

In our recent paper 21 we analyzed conditions that are 
necessary for generation and control of the voltage drop 
across the device. In particular, we focused on the di¬ 
rection and magnitude of transport as well as its de¬ 
pendence on the system parameters. We found the in¬ 
triguing features of a negative absolute and differential 
conductance^^ repeated voltage reversaliPS and noise 
induced voltage reversals. 2 ^ We showed how the direction 
of transport can be controlled by the applied magnetic 
field. Moreover, very recent important aspects con¬ 
cerning the quality of transport occurring in this system 
have been addressed. We analyzed fluctuations of the 
voltage and energetics of the device.^HUl it turned out 
that the power delivered by the external current depends 
not only on its amplitude and frequency but also on ther¬ 
mal noise and the external magnetic field.®®- We explored 
a set of the system parameters to reveal a regime where 
the voltage rectification efficiency is globally maximal. 
We detected the surprising feature of the thermal noise 
enhanced efficiency®^ an d showed how the efficiency of 
the device can be tuned by adjusting the external mag¬ 
netic flux. However, apart from these well investigated 
problems there are still some questions about the trans¬ 
port properties of the system which should be imposed 
and answered. One important example might be a dif¬ 
fusion procesPk 11 ^ of the Josephson phase. Therefore in 
this paper we focus on the connection between the di¬ 
rected transport quantified by the averaged voltage (V), 
its quality measured by the efficiency rj of the device and 
diffusion of the Josephson phase ip. 


III. TRANSPORT QUANTIFIERS 


It is convenient to convert Eq. ([ 2 ]) into its dimension¬ 
less form. There are several forms of such an equation 
depending on time scaling. Here, we follow Ref. |T5] and 
define the new phase x and the dimensionless time t as 


Ifi + 7T 


t = 


h 


( 6 ) 


2 ’ t c e.RJ[ 

In these new variables, Eq. § reads 

Cx{i)+x{i) = — U'(x(i))+F+acos(uii) + V2D£(i), (7) 


where the dot and prime denotes a differentiation with 
respect to the dimensionless time t and the phase x, re¬ 
spectively. The dimensionless capacitance C is the ratio 
between two characteristic time scales C = 77 / 77 , where 
the relaxation time is 77 = RC. Other re-scaled parame¬ 
ters are F = I0/J1, a = A/J± and w = Hr c . The periodic 
potential U(x) of period 27r takes the fornPH 

U(x) = — sin(x) — ^ sin(2:r + $ e — 7r/2), (8) 

where j = J%/ J\. This potential is symmetric if there ex¬ 
ists xq such that U(xq+x) = U(xo~x) for any a;. If j ^ 0, 


it is generally asymmetric and is called a ratchet type, 
see Fig. [2j However, even for j 7 ^ 0 there are certain val¬ 
ues of the external flux <f> e for which it is still symmetric. 
The rescaled zero-mean Gaussian white noise £(f) has the 
auto-correlation function (£(i)£(s)) = S(t — s) and its in¬ 
tensity D = eksT/hJ\ is the quotient of the thermal and 
the Josephson coupling energy. From now on we will use 
only the dimensionless quantities and therefore we skip 
all hats appearing in ([7|. 

There are several quantifiers characterizing transport 
properties of the system. The most important are the 
current-voltage curves in the asymptotic limit of long 
times when all effects due to initial conditions and tran¬ 
sient processes have quiet down. They can be obtained 
from Eq. ([ 5 ]), the dimensionless form of which reads 35 

/>£+27t /oj 

(v) = lim — / E[x(s)]ds, (9) 

t —>00 27T 

where E[x(s)] denotes averaging over initial conditions 
and all realizations of thermal noise. The stationary di¬ 
mensional voltage is then given as 

(V) = RJi(v). (10) 

The long time average voltage (v) represents the basic 
transport measure. However, there are other transport 
quantifiers like the voltage variance 

= (v 2 ) ~ (v) 2 , ( 11 ) 

which describes voltage fluctuations around its average 
value (v). The voltage drop v(t) across the SQUID typi¬ 
cally ranges within the interval of standard deviations, 


v(t) £ [(u) — 07 , (v) + 07 ]. ( 12 ) 


It means that when 07 > | (v) | the instantaneous voltage 
v(t) may assume the opposite sign to the average voltage 
(v) and transport in not effective. 

Next, we can introduce a measure for the efficiency of 
the SQUID in terms of the rectification of thermal fluc¬ 
tuations. It is defined as the ratio between the energetic 
output and the input power ([F + A cos(wf)]u(f)). This 
quantity follows from an energy balance of the underly¬ 
ing Langevin equation fl7[ ) and corresponds to the well 
known Stokes efficiency ® 36 37 


W 2 _ W 2 

(v) 2 + - D/C (v 2 )~D/C' 


( 13 ) 


One should note that this definition yields a non vanish¬ 
ing rectification efficiency even when there is no external 
static current, F = 0. Furthermore, it follows that this 
result is accordant with our intuition: a decrease of the 
voltage variance cr 2 leads to an increase of the rectifica¬ 
tion efficiency. Consequently, to optimize the effective¬ 
ness of the device one should seek for regimes that maxi¬ 
mize the directed voltage and minimize its fluctuations. 
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FIG. 3. Two sets of illustrative trajectories of the phase motion x(t) across the asymmetric SQUID device. Parameters are 
C = 0.496, a — 1.55, ui = 0.406, 4> e = 7t/2, j = 0.5. In panel (a) and (b) the thermal noise intensity is set to D = 10 -5 and 
D = 0.1, respectively. A distinct different diffusive behaviour is observed. 


The next quantifier characterizing transport is related 
to spread of trajectories of the stochastic phase x(t). It 
is the mean square displacement® of the phase defined 
as 

(A x 2 {t)) = ([x{t) - (z(f))] 2 ) = {x 2 (t)) - {x{t)) 2 . 

One can expect that in the long time regime it grows 
according to a power lav® 

(A x 2 (t)) ( ~ 2 D x t a , (14) 

where the constant prefactor D x is sometimes called a 
diffusion coefficient and the exponent a characterizes a 
type of diffusion®® subdiffusion for 0 < a < 1 , normal 
difussion for a = 1 and superdiffusion for a > 1. Another 
special case is ballistic diffusion for a = 2. When a = 1 
the diffusion coefficient can be determined as 

D x = lim D x (t) = lim ( 15 ) 

£—>■00 t—>00 Zt 

Otherwise the above definition is not constructive be¬ 
cause such a quantity is either zero (subdiffusion) or di¬ 
verges to infinity (superdiffusion). The diffusion coeffi¬ 
cient measures the spreading of trajectories x(t) around 
its mean value ( x{t )), see Fig. [3j Intuitively, when it is 
small then the spread of trajectories is small and trans¬ 
port is more optimal. 

Another way of introducing the diffusion coeffi¬ 
cient i s based on the generalized Green-Kubo response 
theory.®®® The asymptotic long time voltage autocor¬ 
relation function 

C(t,r) = (Av{t)Av(t+r)), Av(t) = v(t)-(v(t)), (16) 


The diffusion process is characterized by the low fre¬ 
quency part of the power spectrum of the voltage 
fluctuation^® 

1 r + °° 

X(s) = - dr e isT C(r) ~ s 1_ “. (19) 

" J — OO 

The diffusion coefficient is connected with the above 
equation via the relatior® 


1 r + °° 

Ac=x(0) = 2 drC(r). 


( 20 ) 


Formulas (19) and (20) are particularly important from 


the experimental point of view as they allow for conve¬ 
nient measurement of both the nature of diffusion and if 
needed also its coefficient. 

Finally, the ratio D x /2-k can be considered as a velocity 
describing the normal phase diffusion over one period of 
the potential U(x). Its relation to the average voltage 
(v) determines the dimensionless Peclet number defined 

a j29l3QI32l 


Pe = 


2n\(v)\ 


( 21 ) 


A large Peclet number indicates a motion of mainly regu¬ 
lar nature. If it is small then random or chaotic influences 
dominate the dynamics. 


IV. REGIME OF GLOBALLY MAXIMAL EFFICIENCY: 
NORMAL DIFFUSION 


can be obtained experimentally. Due to presence of the 
periodic driving of frequency ui this function is periodic 
with respect to the first argument® 

C{t, T ) = C(t + T,t) (17) 

where T = 27r/w is a rescaled period of the ac current. 
Therefore we introduce the time average of the autocor¬ 
relation function C(t,r ), namely 

C(t) = ^ J dtC(t,r). (18) 


We have integrated the Langevin equation 0 by em¬ 
ploying a weak version of the stochastic second order pre¬ 
dictor corrector algoritlmi 42 with a time step typically 
set to about 1CT 3 • 2ir/uj. Since Eq. ([T]) is a second- 
order differential equation, we have to specify two initial 
conditions x(0) and i(0). We have chosen phases a:(0) 
and dimensionless voltages i(0) equally distributed over 
interval [0, 27 t] and [2,2], respectively. All quantities of 
interest were ensemble-averaged over 10 3 10 4 different 

trajectories which evolved over 10 3 10 4 periods of the 
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FIG. 4. Impact of the ac current amplitude a on transport characteristics, (a): the dc voltage drop (v) across the SQUID, 
(b): the Stokes efficiency r/, (c): phenomenological power exponent a describing a type of phase diffusion and (d): the diffusion 
coefficient D x in the regime when a « 1. Other parameters correspond to the regime for which the rectification efficiency r/ 
takes its globally maximal value, namely C = 0.496, a = 1.55, c o = 0.406, D = 10~ 5 , $ e = 7t/2, j = 0.5 and F — 0. 


external ac driving. Numerical calculations were done by 
use of a CUDA environment implemented on a modern 
desktop GPU. This scheme allowed for a speed-up of a 
factor of the order 10 3 ti mes a s compared to a common 
present-day CPU method)- 43 * 44 * 

The system described by Eq. 0 has a 7-dimensional 
parameter space {C, a, w, F, j, <J> e , D}. Below, we study 
a nontrivial ratchet effect by putting the dc current F = 
0. We have performed scans of the parameter space at 
a high resolution to determine the general behaviour of 
the system. The conditions that are necessary for the 
generation and control of the direction of transport in the 
SQUID have been extensively studied in these regimes 
in our previous work. 2 3 Moreover, very recently we also 
reported a regime for which the rectification efficiency 77 
is globally maximal ? 9 In this regime the capacitance is 



FIG. 5. The dc voltage bifurcation diagram in the determin¬ 
istic limit D = 0 as a function of the ac driving amplitude a. 
The remaining parameters are given in Fig. [4] 


C ss 0.496, the amplitude of the ac current reads a « 1.55 
and its is frequency w ~ 0.406. Let us now focus on the 
connection between the voltage drop (v ), its rectification 
efficiency 77 and the phase diffusion process across the 
SQUID in this prominent regime. 

In Fig. [4] we present influence of the ac driving am¬ 
plitude a on all characteristics describing quality of the 
transport process observed in the SQUID. In particular, 
panel (a) shows the dc voltage drop (v) across the device. 
There is a finite window of the parameter a for which the 
directed transport of about equal phase velocity (v) is 
observed. When a is smaller than a ~ 1.53 then the rock¬ 
ing mechanism is too weak to induce the non-negligible 
voltage drop. Similar situation can be discovered for the 
amplitudes greater than a « 1.82. These facts have their 
further consequences in the dependence of the Stokes ef¬ 
ficiency 77 on the ac driving amplitude a. Since this quan¬ 
tity is proportional to ( v ) 2 it vanishes when there is no 
directed transport, see panel (b) of the same figure. The 
most intriguing feature of this plot is emergence of the 
rapid maximum for a ~ 1.55. It is caused by the evi¬ 
dent local minimum of the voltage fluctuations a v in the 
vicinity of this point. We refer the reader to Ref. jT51 for 
a detailed study on this effect. Panel (c) depicts the phe¬ 
nomenological power exponent a characterizing the dif¬ 
fusion process as a function of the ac driving amplitude 
a. It was computed from the slopes of time evolution of 
the mean squared displacement (A x 2 (t)). Surprisingly, 
apart from two clearly visible intervals of the ac driving 
amplitude where a ss 1 there is no phase diffusion. Es¬ 
pecially, it is noteworthy that a wide window of a can be 
observed where there is a finite dc voltage drop (v) ^ 0 
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FIG. 6. Impact of variation of the external magnetic flux <5> e on relevant characteristics in the form of: (a) the directed transport 
measured as the dc voltage drop (v) across the SQUID, (b) the Stokes efficiency rj, (c) phenomenological power exponent a 
describing phase diffusion in this setup and finally (d) the normal diffusion coefficient D x . Other parameters are the same as 
in Fig. [4] 


and simultaneously the power exponent vanishes a = 0, 
yielding in a directed and non-diffusive transport across 
the device. In the neighbourhood of the maximal effi¬ 
ciency for a ~ 1.55 the exponent a is guaranteed to be 
unity. In consequence, the phase diffusion is normal, the 
diffusion coefficient D x has well established physical in¬ 
terpretation and can be conveniently computed by use 
of the formula ( |15| ). Its variation is shown in the panel 
(d) of Fig. s The most important finding is that trans¬ 
port in this regime is essentially stochastic. However, the 
diffusion coefficient is small D x « 0.05. Moreover, it is 
seen that one can change the magnitude of D x by ad¬ 
justing the ac current amplitude a. In particular, a small 
change of a is accompanied by multiple increase of the 
diffusion coefficient D x thus leading to the phenomenon 
of diffusion enhancement 

To gain insight into the nature of the phase diffusion 
process we have computed the corresponding determin¬ 
istic (D = 0) dc voltage bifurcation diagram. In Fig. 
[ 5 ] we show the asymptotic long time voltage (v)t aver¬ 
aged over the period T = 2t:/uj of the external ac driving 
for 1024 different initial conditions x(O) and i(O) ran¬ 
domly sampled from the intervals [0, 27 t] and [—2,2], re¬ 
spectively. Therefore all existing attractors are plotted at 
a given value of control parameter a. This figure reveals 
an unexpected simplicity of the phase space dynamics. 
The most important observation is that the system de¬ 
scribed by Eq. ([TJ) is non-chaotic in the analyzed regime 
of the parameters. Moreover, there are regions for which 
only a single period one attractor exists, meaning that 
eventually all initial conditions evolve to it. Thus, each 
trajectory undergoes the same kind of motion resulting 


in the power exponent a = 0 (cf. Fig. [ 4 ];). Moreover, 
there are also intervals where several attractors coexist. 
This, in fact, is enough to observe the phenomenon of de¬ 
terministic diffusion^ 1038 . At sufficiently high temper¬ 
ature, the system will be typically ergodic with thermal 
fluctuations enabling stochastic escape events among co¬ 
existing deterministic separate attractors. In particular, 
transitions between neighbouring periodic solutions give 
rise to diffusive directed transport. 

Let us now briefly discuss an influence of the external 
magnetic flux 4> e on all previously introduced transport 
measures. From the symmetry considerations of Eq. 0 
it follows that for an arbitrary integer number n, the 
transformation <F e —> 2nn — <f> e reverses the sign of the 
average voltage (v) —> —(v). This fact can be directly 
observed in panel (a) of Fig. [6] There are two intervals 
where the average voltage drop assumes non-zero values 
which differ only with the direction of transport. Since 
the external magnetic flux alters the effective potential 
experienced by the phase x we conclude that in order to 
detect the average voltage one must tune it to a given 
rocking mechanism in the form of the ac driving of the 
amplitude a and the frequency u>. Similarly to the pre¬ 
viously discussed case this fact is reflected in the depen¬ 
dence of the Stokes efficiency rj on the external magnetic 
flux $ e . A careful inspection of panel (b) reveals that 
one can tune the effectiveness of voltage rectification just 
by correct adjustment of the external magnetic flux. The 
next plot depicts impact of the external magnetic flux 4> e 
on the exponent a. Regions where the transport is non- 
diffusive dominate this parameter space. It is associated 
with existence of a single period one attractor describing 
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FIG. 7. Influence of variation of the capacitance C on relevant characteristics in the form of: (a) the directed transport measured 
as the dc voltage drop (v) across the SQUID, (b) the voltage fluctuations a v , (c) the Stokes efficiency r/, (d) phenomenological 
power exponent a describing phase diffusion in this setup and finally (e) the normal diffusion coefficient D x . Other parameters 
are the same as in Fig. [4] 


either the running or locked solution of |7| as it was in 
the previous case. Still, there are also intervals where 
the diffusive motion of the phase x can be observed with 
a ss 1. In particular, panel (d) presents the dependence 
of the diffusion coefficient D x on the external magnetic 
flux <3? e for such a scenario: One can conveniently ma¬ 
nipulate the phase diffusion by change of the external 
magnetic flux. This way of control of the diffusion pro¬ 
cess is very convenient from the experimental point of 
view. 

We now focus on the influence of inertia described by 
the capacitance C of the SQUID on the phase diffusion, 
see Fig. [7] Panels (a) and (c) depict the average voltage 
drop ( v) across the device and the Stokes efficiency 77 as a 
function of the capacitance C. They are non-monotonic 
functions similar in shape but without any immediately 
obvious relation to each other. However, the most impor¬ 
tant observation is that in the vicinity of a point corre- 



FIG. 8. The dc voltage bifurcation diagram in the determin¬ 
istic limit D = 0 as a function of the capacitance C of the 
SQUID. The remaining parameters are given in Fig. [4] 


sponding to the maximum transport efficiency rj also the 
voltage drop ( v ) is large. Notably, it lies closely to the 
border between overdamped C —> 0 and damped C ss 1 
regime. In panel (b), fluctuations of the phase veloc¬ 
ity are depicted. One can notice an evident correlation, 
when velocity fluctuations are maximal, the efficiency is 
minimal. The next two panels present the diffusive be¬ 
haviour as a function of the same parameter. In panel 
(d) one can see that the transport is essentially diffusive 
a ss 1 for entire interval of the inertial term variance. The 
last plot depicts the diffusion constant D x versus the ca¬ 
pacitance C of the device. In the overdamped limit the 
transport is rather regular as the diffusion coefficient is 
very small. On the contrary, an increase of inertia is ac¬ 
companied by simultaneous grow of the phase diffusion 
with a sharp increase in the vicinity C ss 0.78. One can 
observe that when the average voltage ( v ) and Stokes effi¬ 
ciency 77 are locally minimal and velocity fluctuations are 
maximal then the diffusion coefficient D x takes its maxi¬ 
mum. Consequently, the transport is highly irregular and 
not optimal. Therefore we validate that quantities char¬ 
acterizing the phase diffusion in the device are somehow 
complementary to those usually used in order to describe 
the quality of transport. They give additional informa¬ 
tion which often corresponds well with the one measured 
by the voltage fluctuations ( pTTj ) or the efficiency (13). It 
is intriguing to find a deeper reason of the sharp increase 
of D x in the vicinity C « 0.78. In Fig. [5J we present the 
dc voltage bifurcation diagram in the deterministic limit 
D = 0 as a function of the capacitance C of the SQUID. 
In the vicinity of C « 0.78, a cascade of bifurcations is 
observed and a transition to chaos takes place. Due to 
this fact, the phase diffusion can be observed even in the 
purely deterministic regime and the diffusion coefficient 

















FIG. 9. Panel (a) the directed transport measured as the dc voltage drop ( v) across the SQUID, (b) the Stokes efficiency r/, 
(c) phenomenological power exponent a describing phase diffusion in this setup and finally (d) the normal diffusion coefficient 
D x presented as a function of thermal noise intensity D. Other parameters are the same as in Fig. [4] 


rapidly increases. 

Finally, we examine the influence of thermal fluctua¬ 
tions on the phase diffusion process. The relevant panels 
comparing all discussed quantities are presented in Fig. 
0 One can see that an increase of the thermal fluctua¬ 
tion intensity D leads to monotonic decrease of both the 
average voltage drop ( v) across the SQUID and its volt¬ 
age rectification efficiency rj. Therefore, we conclude that 
it is a regime for which impact of thermal fluctuations is 
destructive. It is confirmed in the last two panels of this 
figure where the phase diffusion process is studied. Not 
unexpectedly, the system is diffusive for the entire range 
of thermal fluctuation intensity D. It is because thermal 
noise activates stochastic transitions between coexisting 
deterministic disjoint attractors. Moreover, in the last 
panel we can see that when the temperature grows the 
transport becomes more and more diffusive as the diffu¬ 
sion coefficient D x nronotonically increases. 

We collect part of our results by presenting the depen¬ 
dence of the dimensionless Peclet number on the SQUID 
capacitance C and the thermal noise intensity D. The 
corresponding panels can be found in Fig. 10 In the first 
one the reader can observe clear separation between the 
overdamped C —> 0 and underdamped C«1 regimes. In 
the former case the Peclet number is large indicating the 
transport of predominantly regular nature. Contrary, in 
the latter one it is small saying that the phase motion 
is chaotic or diffusive. Having in mind that this partic¬ 
ular regime is essentially non-chaotic, we conclude that 
in that region of a parameter space transport is diffu¬ 
sive. This fact support our previous statement on that 
subject. Moreover, the most important remark is that 
the Peclet number corresponding to the point C = 0.496 


of globally maximal Stokes efficiency rj is relatively large 
being evidence of the diffusive but still highly regular 
transport. In panel (b) of the same figure we can ob¬ 
serve a fast nronotonic decay of the Peclet number with 
increasing temperature of the system. Therefore, in this 
regime thermal fluctuations have destructive impact on 
all transport quantifiers, starting from the average volt¬ 
age ( v ), by the Stokes efficiency ?y of the SQUID and 
finally up to irregularity of the phase motion reflected in 
the large diffusion coefficient D x . 


V. SUMMARY 

We have studied diffusion process of the Josephson 
phase in the asymmetric SQUID system. Our analysis 
has been restricted to the regime of the maximal Stokes 
efficiency of the device and found normal diffusion with 
the exponent a ~ 1. We have searched a neighbourhood 
of this point in the parameter space to check robustness of 
the observed behaviour in domains {a, <f> e , C, D}. When 
the amplitude a of ac current is changed two windows of 
normal diffusion are detected (panel (c) in Fig. 4). In 
turn when the applied magnetic field is varied four inter¬ 
vals of normal phase motion are observed (panel (c) in 
Fig. 6). These areas are interrupted by windows where 
the exponent a = 0 and there is no diffusion. The depen¬ 
dence on the capacitance C and intensity D of thermal 
fluctuations is robust. The exponent a ~ 1 for over¬ 
damped and damped case C < 1 (panel (d) in Fig. 7) 
and rescaled temperature D which can change several 
orders of magnitude (panel (c) in Fig. 8). 

We have presented the possibility of convenient ma- 
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FIG. 10. The dimensionless Peclet number is presented as a function of the capacitance C of the SQUID and the thermal noise 
intensity D in panel (a) and (b), respectively. Other parameters are listed in Fig. [4] 


nipulation of the diffusion coefficient D x by tuning the 
experimentally accessible parameters of the setup like the 
ac driving amplitude a or the external magnetic flux <& e . 
Surprisingly, by doing it one can change its value sev¬ 
eral times. This lead us to the phenomenon of diffusion 
enhancement. 

Last but not least, we have found that the regime of 
maximal Stokes efficiency is essentially non-chaotic. Re¬ 
gions where only a single period one attractor exists dom¬ 
inate the parameter space in the vicinity of this promi¬ 
nent area. Consequently, the phase motion is there non- 
diffusive. However, there are also intervals where due to 
the coexistence of several deterministic separate attrac¬ 
tors the phase diffusion process can be observed. Then, 
sufficiently large thermal fluctuations enable stochastic 
transitions between them resulting in the diffusive di¬ 
rected transport. Overall, thermal noise has destructive 
impact on all presented quantities measuring the qual¬ 
ity of transport in the system, starting from the average 
voltage ( v) and ending on the diffusion coefficient D x . It 
is remarkable that the non-diffusive or regular nature of 
the phase motion across the device is detected mainly in 
the overdamped regime C —> 0. Contrary, when the iner¬ 
tial term is increasing then the transport becomes more 
and more diffusive. 

Finally, an interesting question concerns a pos sibility 
of observing the anomalous phase diffusion 38 39 in this 
setup. This question has been answered in the positive 
for the case of the extremely underdamped (hamiltonian) 
symmetric SQUID device in the deterministic limit.® 
Our case is much more complicated, however, with the 
help of the computational power of modern GPU com¬ 
puters it should still be doable. This additional aspect is 
on our agenda for a potential future research. 

The results described above may be helpful for further 
understanding of nontrivial response of nonlinear dynam¬ 
ics to external driving. They can readily be experimen¬ 
tally verified with an accessible setup consisting of three 
resistively and capacitively shunted Josephson junctions 
formed in an asymmetric SQUID device. Having in mind 
the mechanical interpretation of the studied model our 
research may also have potential applications for particle 
mixing, homogenization, selection or separation tasks.® 
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